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Abstract. The elastic scattering of gqq, gqq and gqq and the thermalization of quark- 
gluon matter are studied. According to Feynman diagrams at the tree level, squared 
^ i' amplitudes for the elastic gqq scattering and the elastic gqq scattering are derived 



in perturbative QCD. Transport equations including the squared amplitudes for the 
QJ ' elastic gqq, gqq and gqq scattering are established. Corresponding to anisotropic gluon 

i-J^ . and quark distributions created in central Au-Au collisions at RHIC, solutions of the 

transport equations show that thermalization time of quark matter can be shortened 
by the elastic gqq, gqq and gqq scattering. 
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O ' 1. Introduction 

o ■ 

I 

Scattering takes place in quark-gluon matter that is created in high-energy heavy- 
ion collisions. If quark-gluon matter has a low number density, two-body scattering 
substantially affects the evolution of quark-gluon matter. The two-body scattering 
. includes the 2-to-2 scattering and the 2-to-3 scattering [Il[2l[3llll|5l|6l[7]. If the number 
density is high, the three-body scattering becomes important. This has been shown by 
the elastic gluon-gluon-gluon scattering in gluon matter which has a number density 
of the order of 19 fm~^ that is reached in central Au-Au collisions at the Relativistic 
Heavy Ion Collider (RHIC) [8]. One effect derived from the elastic gluon-gluon-gluon 
scattering is on the rapid thermalization of gluon matter. The elastic 3-to-3 scattering 
of a heavy quark contributes significantly to the heavy quark momentum degradation in 
the quark-gluon plasma [9]. Since at the Large Hadron Collider (LHC) Pb-Pb collisions 
will be carried out to produce quark-gluon matter which has a higher number density 
than matter reached at RHIC, the elastic 3-to-3 scattering is more involved. Therefore, 
we must study the elastic 3-to-3 scattering and its effects. Due to the complication of 
the study ^0\, in this work we are restricted to the elastic gluon-quark-quark scattering, 
the elastic gluon- quark- ant iquark scattering and the elastic gluon-antiquark-antiquark 
scattering, and apply the scattering to thermalization of quark-gluon matter. 
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We display Feynman diagrams at the tree level for the elastic gluon-quark-quark 
scattering in Section 2 and the elastic gluon-quark-antiquark scattering in Section 3. In 
Section 4, as two examples, we show squared amplitudes for a diagram of the elastic 
gqq scattering and for a diagram of the elastic gqq scattering, respectively. We present 
transport equations that include the elastic scattering of gqq, gqq and gqq. In Section 5 
numerical solutions of the transport equations and relevant discussions are given. The 
last section contains the summary. 

2. Elastic gluon-quark-quEirk scattering 

Some Feynman diagrams for the elastic gluon-quark-quark scattering are shown in Figs. 
1-4. The wiggly lines and solid lines stand for gluons and quarks, respectively. The 
other Feynman diagrams are derived from the diagrams in Figs. 1-4 as follows. 

The six diagrams in Fig. 1 lead to six new diagrams by moving the external gluons 
from the left quark line to the right quark line, i.e. by moving the initial gluon from 
an initial (final) quark to another initial (final) quark and moving the final gluon in the 
same way. If the two final quarks are identical, the exchange of the final quarks in the 
above twelve diagrams leads to twelve more diagrams. We can thus derive 18 diagrams 
from the six diagrams in Fig. 1. 

The four diagrams in Fig. 2 lead to four new diagrams by moving the initial gluon 
from the initial (final) state of a quark to the final (initial) state of the same quark 
and moving the final gluon in the same way. If the two final quarks are identical, the 
exchange of the final quarks in the above eight diagrams leads to eight more diagrams. 
We can thus derive 12 diagrams from the four diagrams in Fig. 2. 

For the elastic scattering of one gluon and two identical quarks, we need to take 
into account 40 diagrams that do not contain any triple-gluon vertex. Each of the 40 
diagrams contains one gluon propagator and four gluon-quark vertices. The 40 diagrams 
form the first class of processes for the elastic gqq scattering. 

Every diagram in Fig. 3 contains one triple-gluon vertex. We can derive 18 
diagrams from the six diagrams. Each of the diagrams D_m and D+m contains one gluon 
propagator between the triple-gluon vertex and the left quark. The diagrams D_m and 
D+M lead to two new diagrams by letting the gluon propagating from the triple-gluon 
vertex to the right quark. The diagrams Dguml, Dgumh, Dgdmh and Dgdml lead to 
four new diagrams by moving the initial or final gluon from the left quark to the right 
quark. If the two final quarks are identical, the exchange of the final quarks in the above 
twelve diagrams generates twelve more diagrams. 

While one triple-gluon vertex is involved in the elastic scattering of one gluon and 
two identical quarks, twenty-four diagrams need to be taken into consideration. Each 
of the 24 diagrams contains two gluon propagators and three gluon-quark vertices. The 
24 diagrams form the second class of processes for the elastic gqq scattering. 

The four diagrams in Fig. 4 are characterized by two triple-gluon vertices or one 
four-gluon vertex, and give rise to four more diagrams by the exchange of the two 
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final quarks if the quarks are identical. The eight diagrams contain three or two gluon 
propagators and two gluon-quark vertices. The 8 diagrams form the third class of 
processes for the elastic gqq scattering. 

We have arrived at the three classes of processes represented by the 72 diagrams 
for the elastic scattering of one gluon and two identical quarks. If the two quarks do not 
possess the same flavor, the processes with the exchange of flnal quarks do not happen 
and thirty-six diagrams are needed. To ensure gauge invariance we include elastic ghost- 
quark-quark scattering for which some Feynman diagrams are shown in Fig. 5. Ghosts 
are indicated by the dashed lines. Each of the diagrams D_]v[fp and D+mfp contains one 
gluon propagator between the ghost line and the left quark line. If the gluon propagator 
is connected to the right quark line, two new diagrams are generated by the diagrams 
D-MFP and D+mfp- In total, we need 7 diagrams for distinguishable quarks and 14 
diagrams for indistinguishable quarks. 

3. Elastic gluon-quark-antiquark scattering 

Since quark- antiquark annihilation may happen in elastic gluon-quark-antiquark 
scattering, the scattering involves more Feynman diagrams than the elastic gluon-quark- 
quark scattering. We only display some diagrams in Figs. 6-12, but the other diagrams 
can de derived as follows. 

In each of the six diagrams in Fig. 6 the initial and final gluons are attached to 
the quark line. While the initial and final gluons are attached to the antiquark line, six 
new diagrams are produced. The four diagrams in Fig. 7 lead to four new diagrams by 
moving an external gluon from the final quark to the initial quark and moving another 
external gluon from the initial (final) antiquark to the final (initial) antiquark. Hence, 
there are twenty diagrams of which any contains no self-coupling of gluons and no 
quark- antiquark annihilation. 

The external gluons in Fig. 8 are attached to one or two of the quark lines. Similarly, 
the initial and final gluons can be attached to one or two of the antiquark lines and six 
new diagrams are thus created. In Fig. 9 the four diagrams lead to four new diagrams 
by moving one external gluon from the final quark to the initial quark and moving 
another external gluon from the initial (final) antiquark to the final (initial) antiquark. 
Therefore, twenty diagrams exist in the case of quark- antiquark annihilation and no 
self-coupling of gluons. 

In total, from Figs. 6-9 we have 10 diagrams of which each contains no triple-gluon 
coupling and no four-gluon coupling. The 40 diagrams form the first class of processes 
for the elastic gqq scattering. 

Each of the diagrams E_m and E+m in Fig. 10 contains one gluon propagator 
between the triple-gluon vertex and the quark line. The two diagrams generate two new 
diagrams while the gluon propagator is connected to the antiquark line. The diagrams 
Eguml, Equmh, Eqdmh and Eqdml generate four new diagrams while the external gluon 
irrelevant to the triple-gluon vertex is moved to the antiquark line. Therefore, twelve 
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diagrams correspond to the case in which every diagram has one triple-gluon vertex and 
quark- ant iquark annihilation and creation do not take place. 

The diagram Eqd (Equ) in Fig. 11 contains one gluon propagator between the 
triple-gluon vertex and the initial (final) quark. The diag rams Eqd and Equ give rise 
to two new diagrams while the gluon propagates between the triple-gluon vertex and 
the initial or final antiquark. In the diagram Eqde the initial quark radiates a gluon 
that creates a quark- antiquark pair. The diagram Eqde generates one new diagram 
by moving the gluon to the initial antiquark. In the diagram Eque the initial quark- 
antiquark pair annihilates into a gluon that is absorbed by the final quark. The diagram 
Eque generates one new diagram by moving the gluon to the final antiquark. Four new 
diagrams are given by the diagrams Equmla, Egumha, Egdmha and Eqdmla by moving 
the external gluon irrelevant to the triple-gluon vertex from a quark to an antiquark. 
Finally, sixteen diagrams correspond to the case in which every diagram has one triple- 
gluon vertex and possesses quark-antiquark annihilation and creation. 

In total, from Figs. 10 and 11 we have 28 diagrams of which each contains one 
triple-gluon coupling. The 28 diagrams form the second class of processes for the elastic 
gqq scattering. 

Any of the diagrams in Fig. 12 has two triple-gluon couplings or one four-gluon 
coupling. The 8 diagrams form the third class of processes for the elastic gqq scattering. 

We have arrived at the three classes of processes shown by the 76 diagrams. If the 
quark-antiquark annihilation may happen, all the 76 diagrams must be considered. If the 
annihilation does not occur, thirty-six diagrams arc needed. To satisfy gauge invariance 
we include elastic ghost-quark-antiquark scattering. Some Feynman diagrams for the 
scattering are plotted in Figs. 13 and 14. In Fig. 13 each of the diagrams E_mfp, 
E+MFP, Eqdfp and Equfp has one gluon propagator between the ghost line and the 
quark line. If the gluon propagator is between the ghost line and the antiquark line, 
four new diagrams are obtained from the diagrams E_mfp, E+mfp, Eqdfp and Equfp- 
In the diagram Eqdefp the gluon radiated from the initial quark breaks into a quark- 
antiquark pair. If the gluon is radiated from the initial antiquark, one new diagram 
is generated from the diagram Eqdefp- In the diagram Equefp the gluon from the 
annihilation of a quark-antiquark pair is absorbed by the final quark. If the gluon is 
absorbed by the final antiquark, a new diagram is generated from the diagram Equefp- 
Finally, we have 6 diagrams derived from the six diagrams in Fig. 13 and together with 
the six diagrams in Fig. 14 we need to consider 18 diagrams for the elastic ghost-quark- 
antiquark scattering. 

4. Transport equations 

We establish transport equations for quark-gluon matter which consists of gluons, quarks 
and antiquarks with up and down fiavors. Denote the gluon distribution function by 
fgi where i labels the ith gluon in scattering. Let the distribution functions for the up 
quark, the down quark, the up antiquark and the down antiquark be fui, fdi, fm and 
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/jj, respectively, and we assume that they are identical, 

fui fdi fui fdi fqit (-'-) 

where i labels the ith quark or antiquark in scattering. With elastic 2-to-2 scattering 
and elastic 3-to-3 scattering, the transport equation for gluons is 

1 f d^p2 d^p3 d^p4 .4.4. , . 
= -2eJ (2^(2^(2^^'^) ^ (Pi +P2 -P3 -P4) 

X { Y I r [/pl/.2(l + /.3)(1 + /.4) - fgsfA^ + + /.2)] 

J^gd-^gd \ ^" | -^gu^gu 
X[/,1/,2(1 + /,3)(1 - /,4) - /,3/,4(l + /,l)(l - /,2)]} 

d^P2 d^pz d^pi d^p5 d^pe 



1 r d 
~2E[J (2^ 



{2T:f2E2 {2T:f2Ez {2T:f2E^ {2T:f2E^ {2t:Y2Eq 
X {2t^Y5\pi + P2 + P3 - P4 - P5 - Pe) I A^sMsfl T 

X [/,i/,2/,3(1 + /.4)(1 + /,5)(1 + /.6) - /.4/,5/.6(l + + /.2)(l + fgz)] 

I 9GgQ_f, . > |2 I I A/f |2 I I A/< _ _ |2 I I A// - - l^^l 

2 •'^'■ggu-^ggu \ "T | •'^l-ggd-^ggd \ "r | J^^ggu-^ggu \ T \ J^i-ggd^ggd I J 

X [/,l/32/,3(l + fgA){l + /,5)(1 " f,,) ' fgifgMl + fgl){l + /^s)!! " /.s)] 
+ fi'gf^ I -^guu^guu P +2 -^gud^gud P + | -Mgdu-^gdu + ^ I -^gdd^gdd P 
+ I -^ffuu-^gwu r + I -^gud-^gud T + I -^gdu-^gdu ^ + I -^gdd-^gdd 1^ 

+ ^ I -^swu-s-suu 1^ ~'~2^^ -^gud-^gud T + I -^gdu^gdu H + ^ I ^gdd-^gdd H 
x[/,l/,2/,3(l + /g4)(l - /,5)(1 - /,6) - /,4/,5/,6(l + " /,2)(1 - /,3)]} , 



and the transport equation for up quarks is 
dt 



1 f d^p2 d^p3 d^p4 , .4^4. . 
^-2eJ {2.f2E2{2.f2E,{2.f2E, ^^''^ ^ +P^-P^- P^) 

X {go I Mug^ug ? [fglfg2{l - f,^){l + fg^) - fMl - /,l)(l + fg2)] 

X [/,l/,2(l - /,3)(1 - /,4) - /,3/,4(l - /,l)(l - /,2)]} 

(i-'^P2 d^ps d^p4 d^p5 d^pQ 



2Ei J (27r)32E2 (27r)32E3 (27r)32E4 (27r)32E5 (27r)32S6 

X (27r)^5^(pi + P2 + P3 - P4 - P5 - Pe) I Y I P 

X [/,l/,2/,3(l - /,4)(1 + fg,){l + /,6) " fg4fgjge{l ' + /.2)(1 + /,3)] 



(2) 
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\ + I Mudg^udg V + \M 

udg I / 

X [fMgsil - /,4)(1 - /,5)(1 + /,6) - f,JMl - /,l)(l - /,2)(1 + /.s)] 
S'qIy^ I -^WMW— I -^uud-^uud I I -^Jidu— >-Ma!« I ) ^ I -^udd^udd \ 

+ I -^uuu^uuu I +77 I -^MMj-i-uud I + I -Mudu^udu I + I ■Mndd^udd I 



2 I - u, Lt u, — r LLUiU, I 1 2 



+ . I J^uuu^uuu I -^uud-^uud I ~^ I -^udu^udu I ) ~l~ ^ I 



udd^udd 



X [/,l/,2/,3(l - /,4)(1 - /,5)(1 - /,6) - /,4/,5/,6(l " /,l)(l " /,2)(1 " /.s)]} , 

(3) 

where vi is the velocity of the massless gluon or up quark; the colour-spin degeneracy 
factors are go = 16 for the gluon and gq = 6 for the quark; pi and p2 {pi, P2 and p^) 
denote the four-momenta of the two (three) initial particles, and p^ and p4 (p4, p^ and 
Pe) of the two (three) final particles in 2-to-2 (3-to-3) scattering; and Ei is the energy 
component of pt. The squared amplitudes of order for the elastic 2-to-2 scattering, 
I Adgg^gg p, | Aigu^gu \^ , ctc, cau bc fouud lu Rcfs. [II1[12]. The squared amplitude 
of order for the elastic ggg scattering, | -Mggg^ggg P, was obtained in the work of 
Ref. [H]. Since the elastic gluon-gluon-quark scattering involves a lot more Feynman 
diagrams than the elastic gqq scattering and the elastic gqq scattering, it will take two 
years to derive the squared amplitude for the elastic gluon-gluon-quark scattering. For 
the time being, we have to give up the elastic scattering of both ggq and ggq, i.e. set 

M-ugg^ugg = ■M.ggu^ggu = M.ggd~^ggd = M.ggu^ggu = -^ggd-i-ggd — 0- 

The squared amplitude for the elastic gluon-antiquark-antiquark scattering equals the 
one for the elastic gqq scattering, for example, 

I A/f -- -- 1-^— I A/f 1^ I A/f - - 1-^— I A/f 1^ 

I ■'^^guu—^guu \ \ ■'^'■guu-^guu \ i \ ■'^^ gud^gud I I ■'^^gud—^gud \ 

The squared amplitude | Aiggg^ggg p is obtained from | Mggg^ggg p by the replacement 
of pi ^ ps and | Mggg^ggg ^ from I Mggg^ggq ^ by thc replaccmeut of pi ps, 
P2 — ^ Pi and ps —> P2- The squared amplitudes for the elastic quark-quark-quark or 
antiquark-antiquark-antiquark scattering and for the elastic quark-quark-antiquark or 
quark-antiquark-antiquark scattering were obtained in the work of Ref. |13] and of 
Ref. [H], respectively. Similar equations for down quarks, up antiquarks and down 
antiquarks can be established. 

The squared amplitude, | Mgqq^gqq ioi g{pi)+q{p2) + q{p-i) g{Pi) + q{p^o) + q{p&) 
or I Mgqq^gqq ^ for g{pi) + q{p2) + q{-pz) ^ 9{Pi) + q{p^) + qi-Pe), is the sum of the 
individually squared amplitudes of the diagrams in the three classes and interference 
terms of different diagrams. The sum of the two interference terms of a diagram in the 
first or third class and a diagram in the second class equals zero. Therefore, there is no 
interference between the first or third class and the second class. Any interference term 
between a diagram in the first class and a diagram in the third class has a very long 
expression and is thus not shown here. Examples of the individually squared amplitudes 
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are the spin- and color-summed squared amplitudes for the diagrams D^lm and E^uu 
which, as shown below, include the average over the spin and color states of the initial 
particles, 

1 1 V- I 12 32g8 1,2 2 

3 ^ ' -^D-LM I = -5-Y^(Sl2Ml6%5 - Sl2Ml6%4 " S12U16U26U35 

spins, colors 

+ S12M16M26M34 - Sl2M?6"35 + S12M16M34 - S12S31M35 

+ S12S31U34 + 2si2S3iM26^i35 " S12S31M26 + 2S12S31M24M35 

+ 2Si2S3iM24%4 " 2Si2S3iM24«26 " 2S12S31M16U34 - 2S12S31M16M24 

+ Si2S3iUig - 2S12S31U15M35 - 2si2S3ini5M34 + 2S12S31M15U26 

+ 2S12S31M15U16 — Si2S3]^M35 — Si2S3;^M34 -|- Si2S3^M26 
2 

+ ■Si2S3]^'Uie -|- S12S23S31M35 -|- S12S23S31M34 — Si2S23S3i'U26 

- S12S23S31U1Q - + Si2Ml6«34 + SI2S31US5 

- S^2'^31'^^34 — SI2S31U2Q — 2s^2'53l'W24 + ■Si2'53l'Wl6 
+ 2S12S31M15 + -5 12'^31 '5l2'^23'531 j 

/[Sl2(si2 + MI6 + U26 - U34 - M35)(-S31 " ^23 - ^34 - ^35)]^, (4) 

1 1 V- I 12 32g8 1 ,2 

2^ I -^E^uu I = —^Y^{-S3lUl5U26U35 + S31M15M26 
spins, colors 

+ S31M15M24M35 - S31M15M24 - ■S3iM^5M26 + S31M15M24 - S31M15M26 
+ S31M15M24 - 2S23^il5Ml6«35 + 2S23?^15Ml6«26 + 2S23^^15Ml6M24 

- S23UI5U35 + S23UI5U26 + 523^15^24 - 2S23M15M16 

- ■523^i?5 - ■523531^^15^35 + 523531^15^*26 + ^23531^^15^24 

- 2S23531U15M16 — 2s2353iU^5 — 52353-|^Mi5 — 5i2523Ml5tt35 

+ 5i2523^il5^i26 + 512523^15^24 - 5i2523M?5 " 5i252353lttl5) 
/[m15(531 + Ml5 - M24 - U26 + M35)(-523 " ^24 - ^26 + ^15)]^, (5) 

where gs is the gauge coupling constant and nine variables are defined as Si2 = {P1+P2Y, 

523 = iP2 +P3)^ 531 = (P3 +Pl)^ Ml5 = (Pl - PbY , ^16 = (Pl - PsY , ^24 = (^2 - Pif , 

■^26 = {P2 — PqYi W34 = {p3 — p^Y and M35 = (^3 — p^Y- Squared amplitudes for 2-to-4 
processes with two gluons and four quarks were obtained from helicity amplitudes in 
Ref . |[15j and can also be derived from Fortran code named CompHEP [16] . One recent 
review on on-shell methods of scattering amplitudes in perturbative QCD has been given 
in Ref. [T7]. Momenta of five partons among the six partons are randomly generated 
and the momentum of the other parton is given by energy-momentum conservation. 
With such a set of six momenta numerical values of the expressions in Eqs. (4) and 
(5) agree with numerical results of the 2-to-4 processes presented in Refs. [T6] by 
reversing the momentum of one final parton. Such agreement is also established for any 
other set of momenta randomly generated. 

Hadronic matter and quark-gluon plasma exist below and above the critical 
temperature of the QCD phase transition, respectively. Due to the medium screening. 
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the gauge coupling constant in quark-gluon plasma is smaller than in hadronic matter 
and at temperatures high enough is so small that perturbative QCD can be applied. 
The gauge coupling constant decreases while the temperature increases [iHl HE]. The 
higher the temperature is, the better any perturbative expansion converges. While the 
temperature is near the critical temperature, the perturbative expansion breaks down; 
but no one has determined the breakdown temperature regime. 



5. Numerical solutions and discussions 



Time dependence of distributions of gluons, quarks and antiquarks is determined by 
the transport equations which have the initial condition at t = 0.2 fm/c generated by 
HIJING |20j for central Au-Au collisions at a/saT/v = 200 GeV and expressed in the form 

m 



with 



and 



f{k±,y,r,z,t) 



A, 



1 



167tR\ 



9{kuy)- 



-(^-ttanhj/)2/2A| 

v/27rAi. 



(6) 



k± cosh y 

(27r)3 



dN 



k± cosh y dyd'^k± ' 

where Ra, y, t, z and r are the gold nucleus radius, transverse momentum, rapidity, 
time, coordinate in the longitudinal direction and radius in the transverse direction, 
respectively. The gluon and the quark have different dN / dyd^kj^. One thousand and 
five hundred gluons within —0.3 < z < 0.3 fm and r < Ra are created from the 
distribution by the rejection method. Two hundred and fifty quarks or antiquarks of 
the up or down flavor are created in the same region. 

Scattering of two partons happens when the two partons have the closest distance 
less than the square root of the ratio of cross section for 2-to-2 scattering to tt. The 
cross section for gg — )• gg is 

g^ 9 \l7s^ + mnls{s + nl) , 2s(s + 2/i2,)2 



99^99 



+ 



Q{s + 2^^lf tiUs + fil) 

and the cross section for gq gq or gq — )• gq is 



+ 2(s + 2^^)ln 



s + t^D 



(7) 



91^91 



91^91 



11 14, 

-, + -(. + 2^1,) In ^^^^ 



+ ?4i±^, (8) 

where s is the square of the total energy of two colliding particles in the center-of- 
momentum system; /iD is the screening mass formulated in Refs. [22l |23[ [2l] and is 
used to regularize propagators. The coupling constant ag = gs/47r = 0.5 is taken in 
finding solutions of the transport equations. The cross section for the elastic qq, qq or 
qq scattering can be found in Ref. [13] where the fraction 8/9 in Eq. (5) should be 
replaced by 4/9. 



2 ^2 
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Scattering of three partons occurs if the three partons are in a sphere of which the 
center is at the center-of-mass of the three partons and of which the radius is [H] 

2 ^ ]_ f d^Pi d^P5 d^Pe 
~ mJ (27r)32E4 (27r)32E5 (27r)32E6 

X (27r)^5^(pi +P2+P3-P4-P5- Pe) I -^3^3 r, (9) 

where m = 1 for gud — )■ gud, gud — )■ gud, guu — )■ guu, gdd — )■ gdd, gud — )■ gud or 
gdu — )■ gdu, and m = 4 for guu — )■ (^mm, (^rfc? — )■ gdd, guu — > (yf-uw or gdd — > (jfrfrf. 

As shown by Eq. (6), the momentum distribution in the longitudinal direction is 
very different from the one in the transverse direction, but gluons, quarks and antiquarks 
at the initial time have similar anisotropy. Gluon distribution functions in all directions 
overlap at the time t = 0.68 fm/c which corresponds to a thermalization time of the 
order of 0.48 fm/c and quark distribution functions at t = 1.56 fm/c which gives a 
thermalization time of about 1.36 fm/c. These distribution functions are plotted in 
Figs. 15 and 16, respectively. Solid curves in Figs. 15 and 16 stand for the Jiittner 
distribution 

with T = 0.5 GeV and Xg = 0.3 for gluon matter at t = 0.68 fm/c, and 

with T = 0.3 GeV and Xg = 0.3 for quark matter at t = 1.56 fm/c, respectively. 

The solutions of the transport equations indicate that gluon matter thermalizes 
rapidly and quark matter thermalizes slowly. How fast thermalization is depends 
on squared amplitudes and distribution functions [10]. Calculations in perturbative 
QCD show that | Mgu^gu ?=\ Mgd-^gd ?=\ Mgu^ga P=| Mgd^gd P, then the term 
QqH Mgu^gu ? + I Mgd-,gd ? + I Mgu^gu ^ + | ^gd^gd P) lu Eq. (2) cquals 
24 I Aigu-^gu P which is near go \ Aiug^ug P in Eq. (3). Therefore, variation of the gluon 
distribution function caused by elastic scattering of both gq and gq is near variation of 
the quark distribution function caused by the same scattering. Numerical calculations 
lead to a similar conclusion about variation of the gluon and quark distribution functions 
caused by elastic scattering of gqq, gqq and gqq. Therefore, the difference between the 
change of gluon distribution and the change of quark distribution or the difference 
between the thermalization time of gluon matter and the one of quark matter is mainly 
given by the elastic ggg scattering and the elastic qqq scattering. Then we need to see 
the point of view from the ggg and qqq scattering. This is accomplished by the four 
aspects: (1) the gluon distribution function fgi is about 2 times the quark distribution 
function fgf, (2) the maximum of | J^gg^gg p fgifgj in Eq. (2) is an order of magnitude 
larger than that of | Muu^uu P + | Mud^ud ? + \ Muu^uu P + | M^d-^^d fgifgj 
in Eq. (3); (3) the maximum of | M.ggg^ggg P is two orders of magnitude larger than 

that of I Mqqq^qqq P, | Mqqq^qqq ^ Or | Mqqq^qqq (4) thc faCtOr g^ fgi fgj fgk / 12 IS OVer 

four times the factor gqfqifqjfqk- Finally, we understand that the rapid thermalization 
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of gluon matter and the slow thermalization of quark matter result mainly from the fact 
that the elastic gg (ggg) scattering has a larger squared amplitude than the elastic qq 
or qq {qqq or qqq) scattering and gluon matter is denser than quark matter. 

The squared amplitude for the elastic gq {gqq or gqq) scattering is comparable to 
the one for the elastic qq or qq {qqq or qqq) scattering. To get a clear understanding of 
contributions of new terms of the elastic scattering of gqq and gqq in Eqs. (2) and (3), 
we approximate the factors, 1 + fgi and 1 — /^j, by 1. In Eq. (2) the maximum of the 
new term gQ{\ Mgu~,gu P + | Mgd^gd P + | Mgu^gu P + | Mg^-^gd ?){fgifq2 - fgsfqi) 
is about half of the maximum of the term ^ | Aigg^gg P {fgifg2 — fgsfg'i), and the 
maximum of the new term g'^[\ \ Mguu^guu ? +\{\ Mgud^gud P + I Mgdu-,gdu ?) + I \ 

■Mgdd^gdd P + I M-guu^guu P + I M-gud-tgud P + I M.gdu^gdu P + I M.gdd^gdd P +i I 
M guu—^guu \ ~l~2(l "^^gud^gud I ~l~ I -^^gdu^gdu n + 4 I J^gdd^gdd \'^]{fglfq2fq3 — fgAfq5fq6) 

2 

is about one-fifth of the maximum of | Mggg^ggg P {fgifg2fg3-fgifg5fg6)- Therefore, 
the new terms provide small contributions to thermalization of gluon matter. In Eq. 
(3) the maximum of the new term go \ J^ug-^ug P {fqifg2 — fqsfgi) is about 2 times 
the maximum of the term 5(q(| | Muu^uu P + | Mud^ud P + I Muu^uu P + I 
M^d^ud ?){fqifq2 - fqzfqi), sud thc maxlmum of the new term gQgG{\ \ Muug^uug P 

+ I -M-udg-yudg P + I M-uug^uug P + | -M-udg^udg \'^){fqlfq2fg3 — fq4fq5fg6) IS near the 

maximum oi ggl^^ I M uud-^uud I ~l~ I -^^udu^udu P) + 4 I -^udd^udd P 

J-L \ \A - _|2_|_i|A/< - -|2_|_|A/f _ _|2_i_|A/f _ -|2_|_l|A/f__ __|2 
'2 I ■'^'■uuu—^uuu I ~r 2 I ■' '■uud^uud I ~r I ■'^'■udu^udu \ ~r | J^'-udd-^udd I ' 4_ \ ■' '■uuu^uuu \ 

+1(1 M^nd^u^d ? + \M 

udd—^udd Wfqlfq2fq3- fqifqbfqd)- HeUCC, the 

new terms have comparable contributions to thermalization of quark matter. Governed 
by the elastic scattering of qq, qq, qqq, qqq and qqq, a thermalization time of the 
order of 1.55 fm/c was obtained in Ref. [l^ for quark matter with the same initial 
distribution as Eq. (6). Hence, the elastic scattering of gq, gqq and gqq shortens 
the thermalization time of quark matter by the amount 0.19 fm/c. About half of the 
amount is a consequence of the elastic gqq and gqq scattering. The elastic gq scattering, 
the elastic gqq scattering and the elastic gqq scattering have different contributions in 
shortening the thermalization time. Taking up-quark matter as an example, amounts 
by which the thermalization time is shortened by relevant elastic scattering are listed in 
Table 1. The amount by which the thermalization time is shortened by the elastic gqq 
scattering is larger than by gqq. 

In Ref. [8] one thousand gluons are generated from a distribution that is 
homogeneous in coordinate space but anisotropic in momentum space. Gluon matter 
is controlled to evolve in the longitudinal direction and is governed by the elastic ggg 
scattering. The resultant fugacity and temperature are 0.065 and 0.75 GeV, respectively. 
In Ref. [13] six hundred and sixty-six quarks are generated from a distribution similar to 
that in Ref. [8]. Quark matter is also controlled to evolve in the longitudinal direction 
and is governed by the elastic qqq scattering. The resultant fugacity and temperature 
are 0.04 and 0.59 GeV, respectively. Due to the restriction of longitudinal expansion, 
the two fugacities are small and the two temperatures are high. In Ref. [1^ five hundred 
quarks and five hundred antiquarks are created from the same distribution as that in 
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Eq. (6). Governed by the elastic scattering of qqq, qqq, qqq and qqq, quark matter and 
antiquark matter evolve in both the transverse direction and the longitudinal direction. 
The resultant fugacity and temperature are 0.31 and 0.27 GeV, respectively. Without 
the restriction of longitudinal expansion, the fugacities in Ref. [H] and in the present 
work are not small and the temperatures are not higher than those in Refs. [8] and 
[13] . The fugacity and temperature obtained in the present work must be different from 
those in Ref. [H] since the present work involves gluon matter that is absent in Ref. 

Even though the study and application of elastic 4-to-4 scattering is not the 
purpose of the present work, we still can know the occurrence probability of the 
4-parton scattering and it has been shown in Ref. [H]. One thousand and five 
hundred partons were generated from the same distribution as that in Eq. (6) within 
—0.3 < z < 0.3 fm and r < Ra- Given an interaction range of 0.62 fm, the 2-parton 
scattering has the occurrence probability of 30%, the 3-parton scattering 20%, and the 
4-parton scattering 14.6%. Therefore, the elastic 4-parton scattering is expected to give 
a smaller contribution to thermalization than the elastic 3-parton scattering. Interaction 
of partons in a sphere with a radius the same as the interaction range takes place and 
the number of partons in such a sphere is at most 14. This means that at most 14- 
parton scattering is allowed. As a consequence, the occurrence probability of 15-parton 
scattering is zero and the one of the 14-parton scattering is very small. The occurrence 
probabilities for 5-parton, 6-parton, 7-parton, 8-parton, 9-parton, 10-parton and 11- 
parton scattering are 11%, 9%, 7.5%, 4.4%, 2.3%, 0.9% and 0.2%, respectively. The 
occurrence probabilities for 12-parton scattering and 13-parton scattering are negligible. 
The sum of the occurrence probabilities from the 2-parton scattering through the 14- 
parton scattering equals 1 and the occurrence probabilities form a convergent series. 

6. Summary 

We have established the transport equations that include the squared amplitudes for 
the elastic scattering of gqq, gqq and gqq. The squared amplitudes are derived at the 
tree level of the scattering in perturbative QCD and are expressed in terms of the 
nine Lorentz-invariant momentum variables. The elastic scattering of gqq, gqq and gqq 
shortens the thermalization time of quark matter as well as antiquark matter. This 
is an effect of the three-body scattering while the number density is high. In not only 
quark-gluon matter with a high number density but also a very dense scalar field system, 
the elastic 3-to-3 scattering plays a significant role [25] . 
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Figure 1. Elastic gluon-quark-quark scattering. 
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Figure 3. Elastic gluon-quark-quark scattering. 
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Figure 4. Elastic gluon-quark-quark scattering. 
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Figure 5. Elastic ghost-quark-quark scattering. 
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Figure 6. Elastic gluon-quark-antiquark scattering. 
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Figure 7. Elastic gluon-quark-antiquark scattering. 
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Figure 8. Elastic gluon-quark-antiquark scattering. 






Figure 9. Elastic gluon-quark-antiquark scattering. 
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Figure 10. Elastic gluon-quark-antiquark scattering. 
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Figure 11. Elastic gluon-quark-antiquark scattering. 
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Figure 12. Elastic gluon-quark-antiquark scattering. 



Thermalization of quark- gluon matter with the elastic scattering of gqq, gqq and gqq 26 




Figure 13. Elastic ghost-quark-antiquark scattering. 
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Figure 15. Gluon distribution functions versus momentum in different directions 

while ghion matter arrives at thermal equilibrium. The dotted, dashed and dot- 
dashed curves correspond to the angles relative to one incoming beam direction 
6 = 0°,45°,90°, respectively. The solid curve represents the thermal distribution 
function. 




Figure 16. The same as Fig. 15, except for quark distribution functions while quark 
matter arrives at thermal equilibrium. 



